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Abstract 

The most recent data on the weighted transverse single spin asymmetry 
A& m(4> h -<t>s) from HERMES and COMPASS collaborations are analysed within 
LO parton model with unintegrated parton distribution and fragmentation 
functions; all transverse motions are taken into account, with exact kinemat- 
ics, in the elementary interactions. The overall quality of the data is such that, 
for the first time, a rather well constrained extraction of the Sivers function 
for u and d quarks is possible and is performed. Comparisons with models 
are made. Based on the extracted Sivers functions, predictions for A^^^ h ^ 
asymmetries at JLab are given; suggestions for further measurements at COM- 
PASS, with a transversely polarized hydrogen target and selecting favourable 
kinematical ranges, are discussed. Predictions are also presented for Single 
Spin Asymmetries (SSA) in Drell-Yan processes at RHIC and GSI. 



1. Introduction 

In a recent paper pQ we have discussed the role of intrinsic motions in inclusive 
and Semi-Inclusive Deep Inelastic Scattering (SIDIS) processes, both in unpolarized 
and polarized i p — ► I h X reactions. The LO QCD parton model computations have 
been compared with the experimental dependence of the unpolarized cross section 
on the azimuthal angle, around the virtual photon direction, between the leptonic 
and the hadronic planes (Cahn effect [2J); at small transverse momentum of the 
produced hadron h, such an effect is dominantly related to intrinsic motions and it 
allows an estimate of the average values of the transverse momenta of quarks inside 
a proton, k±, and of final hadrons inside the fragmenting quark jet, p ± , with the 
best fit results: 

(k 2 ± ) = 0.25 (GeV/c) 2 (pi) = 0.20 (GeV/c) 2 . (1) 

More detail, both about the kinematical configurations and conventions |3] and 
the fitting procedure can be found in Ref. [Tj. We only notice here that the above 
values have been derived from sets of data collected at different energies and in 
different ranges of the kinematical variables x B , Q 2 and Zh, looking at the combined 
production of all charged hadrons in SIDIS processes; constant and flavour indepen- 
dent values of (kj_) and (p\) have been assumed, avoiding at this stage complications 
related to possible x, z and Q 2 dependences. Rather than a definite derivation, the 
above results are better considered as a consistent simple estimate and a convenient 
parameterization of the true intrinsic motion of quarks in nucleons and of hadrons 
in jets, supported by the available experimental information. 

Equipped with such estimates, in Ref. [T] we have studied the transverse single 
spin asymmetries A s ™^ ^ observed by HERMES collaboration [I]; that allowed 
a first rough extraction of the Sivers function |Hj 

A% p t (x, k ± ) = -^ f^(x, k ± ) , (2) 
nip 

defined by 

f q/p i(x, k±) = f q/p {x, k±) + ^ A N f q/p i(x, S ■ (P x fe ± ) (3) 

= f q/p (x,k ± )-f^(x,k ± ) 5 -( Px ^) ; (4) 

Trip 

where f g / p (x, k±) is the unpolarized x and k± dependent parton distribution (k± = 
\k±\); m p , P and S are respectively the proton mass, momentum and transverse 
polarization vectors (P and k± denote unit vectors). 

The Sivers function extracted from HERMES data in pQ was shown to be consis- 
tent with preliminary COMPASS data on A s Jr^ h ~^ s ' > obtained on a deuteron target, 
in a different kinematical region jHj. 
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While the preliminary HERMES data 4] offered a definite indication of a non 
zero Sivers effect, their amount and quality were not yet such that an accurate 
extraction of the Sivers functions was possible; that reflects in the values of the 
parameters of the Sivers functions given in Table I of Ref. |]Q, which have large 
uncertainties. 

New HERMES data are now available jjj; they are consistent with the previous 
ones, with much smaller errors. Similarly, COMPASS collaboration published jH] 
their preliminary results. We consider here these whole new sets of HERMES and 
COMPASS data and perform a novel fit of the Sivers functions. It turns out that 
the new data constrain much better the parameters, thus offering the first direct 
significant estimate of the Sivers functions - for u and d quarks - active in SIDIS 
processes. The sea quark contributions are found to be negligible, at least in the 
kinematical region of the available data. 

The modeled and extracted Sivers functions A N f u / p t (x, k±) and A N f d / p t (x, k±) 

are used to compute, and thus predict, the values of Ay^ h ^ expected at COM- 
PASS, for scattering off a polarized hydrogen (rather than deuteron) target, which 
avoids cancellations between the opposite u and d contributions. Suggestions for 
the selection of favourable kinematic regions, where the asymmetry is sizeable, are 
discussed. Similar predictions, with strongly encouraging results, are given for po- 
larized SIDIS processes at JLab. 

Finally, we exploit the QCD prediction jHj 



and compute a single spin asymmetry, which can only originate from the Sivers 
mechanism ^^j, for Drell-Yan processes at RHIC and GSI. 

2. Extracting the Sivers functions 

Following Ref. pQ, the inclusive (£p —>■ £X) unpolarized DIS cross section in 
non collinear LO parton model is given by 



DIS — — 



AtO, k ±)\ 



D-Y 



(5) 





and the semi-inclusive one (£p — > £hX) by 




dQ 2 



J-D%{z, P± ), 



(7) 



where 





and 



da £q 
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(9) 



dQ 2 



e " s 2 Q 4 



2 



Q 2 , x B and y = Q 2 /(x B s) are the usual leptonic DIS variables and Zh,Pr the 
usual hadronic SIDIS ones, in the 7*-p cm. frame; x and z are light-cone momentum 
fractions, with (see Ref. pQ for exact relationships and further detail): 



x = x B + O 



'A-1 



z = z h + 



\Q 2 ) \Q 2 ) 

The elementary Mandelstam variables are given by 



p ± = P T -z h k± + 



'kY 



(10) 



s 2 



u 



y 2 
y 2 



1-4 



y cos if 
k± cos (p 



Q V^y 



♦»») 

hi 



(ii) 

(12) 

where tp is the azimuthal angle of the quark transverse momentum, k± = kj_(cosip, 
sin <p,0). Regarding angle definitions and notations we adopt throughout the paper 
the so-called "Trento conventions" |3] (see also Fig. 3 of Ref. [T]). 

The sin — <ps) weighted transverse single spin asymmetry, measured by HER- 
MES and COMPASS, which singles out the contribution of the Sivers function (J2J), 
is given by: 



Asm{4> h -<ps) 
I\ UT 



(13) 



E 



b s dcj) h d 2 k ± A N f q/ ^ (x, k±) sin(y? - (j>s) 



ur , J — Dg(z,p±) sin((j)h — (j>s) 

CLL^ Z}i 



Jd(p s d(p h d 2 kj_ fq/ p (x, k ± ) 



J-D h (z, P± ) 



We shall use Eq. (|13|) . in which we insert a parameterization for the Sivers functions, 
to fit the experimental data. 

We adopt the usual (and convenient) gaussian factorization for the distribution 
and fragmentation functions: 



and 



f q / p (x,k ± ) = f q (x) 



D h (z,p ± )=Dhz) 



1 



7T<*i> 
1 



-k\/(kl) 



-p 2 J(pD 



(14) 



(15) 



with the values of {k±) and (pj_) of Eq. (£[]). Isospin and charge-conjugation relations 
imply 



Df(z) = Df(z) = Df(z) = Df(z) = D iaY (z) 



Dl{z) = Df(z) = Df(z) = Df(z) = Dntoiz) . 



(16) 
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The integrated parton distribution and fragmentation functions f q (x) and D^(z) 
are taken from the literature, at the appropriate Q 2 values of the experimental data 

D3Q3. 

We parameterize, for each light quark flavour q = u,d, u, d, the Sivers function 
in the following factorized form: 

A% p t(x, k ± ) = 2M q {x) h(k ± ) f q/p (x, k ± ) , (17) 

where 

M q {x) = N q x^(l - xf« , (18) 

a q b q 

h(k±) = WTm! ' (19) 

N q , a q , b q and M (GeV/c) are free parameters. f q / p {x,k±) is the unpolarized 
distribution function defined in Eq. ()14j) . Since h(k±) < 1 and since we allow the 
constant parameter N q to vary only inside the range [—1, 1] so that |A/" g (x)| < 1 for 
any x, the positivity bound for the Sivers function is automatically fulfilled: 

|A% P ,(*,*JI S1 



2f q / P (x, k±) 

We have first attempted a fit of the HERMES and COMPASS data, taking into 
account 4 Sivers functions (for u,d,u and d quarks), for a total of 13 parameters, 
like in Ref. pQ. However, it turns out that the available data are almost insensitive 
to the sea quark (and, in general, small x) contributions, which leads to largely 
undetermined parameters of the corresponding Sivers functions. Indeed, we have 
explicitely checked that various choices of A N f a / p t and A N f^ p t do not significantly 

affect the computation of A s Jr^ h Eq. (|T3|) . in the kinematical regions of the 
performed experiments. We have then neglected the contributions of these functions 
and considered only the contributions of A N f u / p T and A N f d / p T, for a total of 7 free 
parameters: 

N u a u b u N d a d b d M . (21) 

The results of our fits are shown in Figs. ^andEI The weighted SSA A s J^^ h 
is plotted as a function of one variable at a time, either Zh or x B or Pt; an inte- 
gration over the other variables has been performed consistently with the cuts of 
the corresponding experiment (see Ref. |]Q for further detail). The resulting best 
fit values of the parameters are reported in Table I. The shaded area in Figs. 
and 121 corresponds to one-sigma deviation at 90% CL and was calculated using the 
errors (Table I) and the correlation matrix generated by MINUIT, minimizing and 
maximizing the function under consideration, in a 7-dimensional parameter space 
hyper-volume corresponding to one-sigma deviation. 
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N u = 


0.32 ±0.11 


N d = 


-1.00 ±0.12 


a u = 


0.29 ±0.35 


ad = 


1.16 ±0.47 


K = 


0.53 ±3.58 


b d = 


3.77 ±2.59 


M'i = 


0.32 ±0.25 (GeV/c) 2 


X 2 /d.o.f. = 


1.06 



Table I: Best fit values of the parameters of the Sivers functions. Notice that the errors 
generated by MINUIT are strongly correlated, and should not be taken at face value; the 
significant fluctuations in our results are shown by the shaded areas in Figs. ^ and [3 

In Fig. 1 we also show predictions, obtained using the extracted Sivers functions, 
for 7T° and K production; data on these asymmetries might be available soon from 
HERMES collaboration. 

3. Comparison with models and predictions for SIDIS processes 

The extracted Sivers functions for u and d quarks are shown in Fig. El where we 
plot, for comparison with other results, the first k± moment 

= l* k ^ AA W(*, k ± ) = -At (% (z) • (22) 

The solid line corresponds to the central values in Table I and the shaded area 
corresponds to varying the parameters within the shaded areas in Figs. ^ and 121 
The other curves show results from models or fits to different data [T3*l rjJJ rjj)|, as 
discussed below. 

• The ^-dependences of both A N f u / p i and A N f d / p i - as modeled in Eqs. ([T71 - 
frn|) and shown in Fig. El appear to be rather well determined, keeping in 
mind that the data are essentially confined in the region 0.01 < x B < 0.2. 
We notice that a 13-parameter fit - including u and d contributions - would 
lead to similar results; however, the values of A^/^/pt and A N f^j p ^ within 
their shaded areas, would be consistent with zero in the kinematical region of 
HERMES and COMPASS experiments. That is why we have not considered 
these contributions here. 

• The large-x behaviour of the Sivers functions cannot be fixed by the existing 
data. According to the counting rules of Ref . [12] , and keeping in mind that 
A N f q / p i- originates from the interference between distribution amplitudes with 
different proton helicities [T71 [TH] , one expects the large-x behaviours 

A% pT ~A% pT ~(l-x) 4 . (23) 

JLab data will cover the appropriate region to help checking this prediction. 
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• The dot-dashed line in Fig. E] shows fit I of Ref. [T3], where the qx/mp 
weighted SIDIS asymmetries were fitted and the large N c relation [THl 120] was 
adopted: 

f#{x,k ± ) = -ftf(x,k ± ). (24) 

Notice that their results are in qualitative agreement with ours and that Eq. 
(J24j) naturally turns out to be approximately true in our fit. 

• The dashed line in Fig. El plots the first k± moment of the Sivers functions 
obtained in Ref. [Hj, by fitting An data in p" p — > n X processes; these data 
are mainly sensitive to large x value, where - again - approximately opposite 
values of A N f u / p r and A N f d / p t seem to be favoured. However, as discussed in 
Ref. [T] , the universality of the Sivers functions active in SIDIS and p 1 ' p — > tt X 
processes is still an open issue. 

• Most theoretical models give a Sivers function for u quarks much larger, in 
magnitude, than for d quarks; this can be seen, for example, from the dotted 
curve in Fig. El taken from the computation, in the MIT bag model, of Ref. 
[To] . The same is true for the Sivers functions obtained, within a spectator 
model with diquarks, in Refs. [2T], |22j and [213] . 



3.1 AuY h 9s> at COMPASS with polarized hydrogen target 

By inspection of Eq. ([THJ) it is easy to understand our numerical results for the u 
and d Sivers functions. In fact one can see that for scattering off a hydrogen target 
(HERMES), one has 

(A$*-**>) hydrogen ~ 4 A% pT D h u + A% pT D\ , (25) 

while, for a scattering off a deuterium target (COMPASS), 

«*-«) deuterium ~ (A~/„ /p , + A% pt ) (4D* + D'i) . (26) 

Opposite u and d Sivers contributions suppress COMPASS asymmetries for any 
hadron h. These opposite contributions do not affect the 7r + asymmetry measured 
off a a hydrogen target, Eq. (|2*3|l : in this case the charge factor 4 and the favourite 
fragmentation function (D^ + > D% + ) combine to make the first term of Eq. (J25j) 
larger than the second one. The cancellation between the two terms is stronger 
for the 7T~ asymmetry, because in this case the charge factor 4 in the first term of 
Eq. ()25|) couples to the unfavoured fragmentation function (D^ < D% ). Similar 
arguments hold for the production of kaons and, in general, for the production of 
charged hadrons, which is dominated by pions. 

However, the COMPASS collaboration will soon be taking data with a trans- 
versely polarized hydrogen target. We can easily compute the expected results: 
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adopting the same experimental cuts which were used for the deuterium target pQ 
we obtain the predictions shown in the upper panel of Fig. |3J The asymmetry is 
found to be around 5%. These expected values can be further increased by properly 
selecting the experimental data, thus excluding kinematical regions whose contribu- 
tion to the asymmetry is negligible. For example, selecting events with 

0.4 < z h < 1 0.2 < P T < 1 GeV/c 0.02 < x B < 1 , (27) 

yields the predictions shown in the lower panel of Fig. 0] The asymmetry for 
positively charged hadrons becomes larger, and, provided that enough statistics can 
be gathered, one expects a clear observation of a sizeable azimuthal asymmetry also 
for the COMPASS experiment. 

3.2 Au^ h ^ at JLab with polarized hydrogen target 

Also JLab experiments are supposed to measure the SIDIS azimuthal asymmetry 
for the production of pions on a transversely polarised hydrogen target, at incident 
beam energies of 6 and 12 GeV. The kinematical region of this experiment is very 
interesting, as it will supply information on the behaviour of the Sivers functions in 
the large-x s domain, up to x B ~ 0.6. The experimental acceptance for JLab events 
at 6 GeV is constrained by |24j : 

0.4 < < 0.7 0.02 < P T < 1 GeV/c 0.1 < x B < 0.6 (28) 
0.4 < y < 0.85 Q 2 > 1 (GeV/c) 2 W 2 > 4 GeV 2 1 < E h < 4 GeV. 

while, with an incident beam energy of 12 GeV, this becomes: 

0.4 < z h < 0.7 0.02 <P T < 1.4 GeV/c 0.05 < x B < 0.7 (29) 
0.2 < y < 0.85 Q 2 > 1 (GeV/c) 2 W 2 > 4 GeV 2 1 < E h < 7 GeV. 

Imposing these experimental cuts we obtain the predictions shown in Fig. El 
A large and healthy azimuthal asimmetry for 7r + production should be observed. 
Similar results have been obtained also in an approach based on a Monte Carlo 
event generator [23]. However, one relevant comment is in order: 

• As the region of high x B is not covered by HERMES and COMPASS experi- 
ments, the predictions for the large-x B dependence of the asymmetry are very 
sensitive to the few large-x s data points of these two experiments. As a con- 
sequence, the results for JLab experiments may still change drastically in the 
region 0.4 < x B < 0.6, and the asymmetry might be much smaller than pre- 
sented in Fig. This reflects in the wide shaded area at large x B values. 
Conversely, the results on Py and Zh dependences are more stable as they de- 
pend on the x s -integrated Sivers function. Notice also the little dependence 
on the beam energy, consistent with the approximate factorized form of the 
numerator and denominator of Eq. (fTH|) , which leads to cancellations in their 
ratio. 
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4. Transverse single spin asymmetries in Drell-Yan processes 

Let us now consider the transverse single spin asymmetry, 

da^ — da^ . . 

A " = ■ < 3 °> 

for Drell-Yan processes, p 1 ' p — > £ + £~ X, p^ p —>■ £ + £~ X and p 1 * p — > £ + £" X, where 
da stands for 

d 4 a 



dydM 2 d 2 q^ ^ 

and y, M 2 and q T are respectively the rapidity, the squared invariant mass and the 
transverse momentum of the lepton pair in the initial nucleon cm. system. The 
cross section can eventually be integrated over some of these variables, according to 
the kinematical configurations of the experiments. 

In such a case the single spin asymmetry (fHU|) can only originate from the Sivers 
function and is given (selecting the region with q\ «C M 2 , q? — k±) by 



An = E q e 2 , J d 2 k ±q d 2 k ± g 5 2 (k ±q + k ±q - q T ) A N f q/ ^ (x q , k ±q ) fg/ p (x q , k±g) 
2J2 q e 2 f d 2 k_L q d 2 k ±q S 2 (k ±q + k ±q - q T ) f q / P (x q , k ±q ) fg/ P (x q , k ±q ) 

where q — u,u, d, d, s, s and 

y/S y/S 

Eq. (jH2| explicitely refers to p 1 ' p processes, with obvious modifications for p^ p and 
p^ p ones. 

Inserting into Eq. (|32|) the Sivers functions extracted from our fit to SIDIS data 
and reversed in sign according to Eq. ©, we obtain the predictions shown in Figs. 
6 and 7. Fig. 6 shows the value of A N as a function of M and xp = x q — x q , for 
RHIC configurations: the lepton pair transverse momentum q T has been integrated 
in the range < q^ < 1 GeV/c, while the rapidity variable y and the lepton pair 
invariant mass M have been integrated according to the experimental situations, as 
indicated in the legenda. The integration over the azimuthal angle of q T has been 
performed, to avoid cancellations, as in Ref. ^Oj, that is integrating over <pq T in the 
range [0, 7r/2] only, or, alternatively, taking into account the change of sign in the 
different production quadrants. In either case, the integration gives an overall 
factor 2/n. As the shaded areas in previous figures, the closed areas correspond to 
the uncertainty in our determination of the Sivers functions. 

Fig. 7 shows the same plots for the PAX experiment planned at the proposed 
asymmetric pp collider at GSI: q T has been integrated over the same range as for 
RHIC predictions, while y and M as indicated in the legenda. Results for p^ p and 
p^ p processes are identical, due to charge conjugation invariance. Notice that in 
our configuration the polarized proton or antiproton always moves along the {+z)~ 
direction. 

The correct interpretation of these results requires some further considerations. 
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• In our computations we have used the value of (k\) = 0.25 (GeV/c) 2 , obtained 
from an analysis of SIDIS data pQ; such a value is certainly appropriate for 
consistently computing spin asymmetries in SIDIS processes, in the 7*p cm. 
frame, as we have done in Section 3.1 and 3.2. This value naturally corresponds 
to the intrinsic motion of partons confined in a nucleon, simply according to 
uncertainty principle, and describes well the P? dependences of measured 
cross sections, up to Pt — 1 GeV/c. In addition, as we have seen, it allows an 
understanding of the azimuthal asymmetries, which would otherwise vanish. 

• However, when considering other processes, as the inclusive production of 
hadrons or leptons in pp or pp interactions, we know that higher order QCD 
corrections, like the threshold resummation of large logarithms due to soft 
gluon emission |2Zj , lead to large .fT-factor enhancements of the cross sections. 
Our k± unintegrated approach to the description of hard scattering processes 
within a generalization of the QCD factorization theorem |T3J EH]- can be 
considered as an effective model which not only takes into account the original 
partonic intrinsic motion (related to parton confinement), but also, to some 
extent, the intrinsic k± built via soft gluon emission. Indeed, the values of 
(kj_) used in Ref. (Hj in order to describe the data on the unpolarized pp — ► 
7r X processes are higher than the values used here, and those requested for 
the Drell-Yan cross-section might be even higher. The average (k±) estimate 
of 0.25 (GeV/c) 2 might be at most adequate to explain the Drell-Yan cross 
section up to qx < 1 GeV/c, but would badly fail above that value. 

• For the reasons explained above, consistently with our approach expected to 
hold in the k± ~ Pt — qr region, in our predictions for A^, Figs. 6 and 7, 
we have integrated over q T up to qx — 1 GeV/c. In addition, we notice that 
the value of A N , as given by Eq. (|32|h is little sensitive to the chosen value of 

while both the numerator and the denominator of Eq. ()32|) greatly vary 
with (k\), their ratio does much less so. We then consider our predictions for 
An, assuming the validity of the relation (jSJ), safe and significant. 

5. Comments and conclusions 

We have considered the most recent data from polarized SIDIS processes which 
single out the Sivers effect, namely the A s ^^ h ^ transverse single spin azimuthal 
asymmetry, measured by HERMES [7J and COMPASS jH] collaborations for charged 
hadron production. Assuming a Gaussian factorization of the k± and p± dependence 
of all distribution and fragmentation functions, together with a most simple param- 
eterization of the ^-dependence of the unknown Sivers functions, we have exploited 
the data to extract information on A N f u / p ^(x, k±) and A N f d / p t(x, k±). 

For the first time, the amount and quality of the experimental results allow a 
significant, although still limited in x-range, estimate of the Sivers functions for u 
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and d quarks; these turn out to be definitely different from zero, well inside the 
positivity bound of Eq. (j2*Uj) and almost opposite to each other. This last feature, 
predicted theoretically in some models [131201, explains naturally and is related to 
the small asymmetry observed by COMPASS in scatterings of muons off a deuteron 
target. 

According to the general strategy of combining new experimental information 
with the computation and prediction of new expected results, the extracted func- 
tions have been used to compute A s ^^ h ^ in other experiments. It turns out 
that, contrary to the results obtained off a polarized deuteron target, a sizeable h + 
asymmetry should be measured by COMPASS collaboration once they switch, as 
planned, to a transversely polarized hydrogen target; a careful choice of the kine- 
matical region of the selected events would help in further increasing the numerical 
value of the asymmetry for positively charged hadron production. 

Large values of A s ^p^ h ~^ s ' are expected at JLab, both in the 6 and 12 GeV oper- 
ational modes, for 7r + inclusive production; in particular the Zh and Pt dependence 
of the asymmetry seems to be stable and reliable, while the x B dependence shows 
large uncertainties due to the lack of HERMES and COMPASS information in this 
region. JLab experiments have the unique features of exploring the large x B be- 
haviour of the quark distribution functions, where predictions from QCD counting 
rules, Eq. (J2BI, could be tested. 

The Sivers effect was believed for some time to be forbidden by QCD time 
reversal properties [2H|; however, this proved to be incorrect jS] after an explicit 
model showed the existence of a non zero Sivers function |18j . The original proof 
of the vanishing of the Sivers effect turned into the relation of Eq. (jSJ), which 
predicts opposite values for the Sivers functions measured in SIDIS and Drell-Yan 
processes. We have then used this basic QCD relation and computed the single spin 
asymmetries in Drell-Yan processes given in Eq. (|32|) : these can only be generated 
by the Sivers functions, since no fragmentation functions are needed to describe 
this process. We have used the same functions as extracted from SIDIS data, with 
opposite signs. The predicted could be measured at RHIC in pp collisions and, 
in the long range, at the proposed PAX experiment at GSI [26j, in pp interactions. 
It would provide a clear and stringent test of basic QCD properties. 

A phenomenological study of the Sivers asymmetry - the correlation between 
the intrinsic k± of partons and the proton spin - is now possible, thanks to the 
existing experimental information and more which will soon be available. Basic 
properties of the QCD proton structure can and will be clarified. A good control 
of the Sivers mechanism will help in learning and understanding about other fun- 
damental partonic spin properties, like the transversity distribution [3U] and the 
Collins mechanism 129 . 



Acknowledgements 

We would like to thank Elke Aschenauer, Harut Avakian and Delia Hasch for 




10 



fruitful discussions. We acknowledge the support of the European Community- 
Research Infrasctructure Activity under the FP6 "Structuring the European Re- 
search Area" programme (HadronPhysics, contract number RII3-CT-2004-506078). 
U.D. and F.M. acknowledge partial support by MIUR (Ministero dell'Istruzione, 
dell'Universita e della Ricerca) under Cofinanziamento PRIN 2003. 



References 

[1] M. Anselmino, M. Boglione, U. D'Alesio, A. Kotzinian, F. Murgia and A. Prokudin, 
Phys. Rev. D71 (2005) 074006 

[2] R.N. Cahn, Phys. Lett. B78 (1978) 269; Phys. Rev. D40 (1989) 3107 

[3] A. Bacchetta, U. D'Alesio, M. Diehl and C.A. Miller, Phys. Rev. D70 (2004) 117504 

[4] HERMES Collaboration, A. Airapetian et al, Phys. Rev. Lett. 94 (2005) 012002; 
U. Elschenbroich, G. Schnell and R. Seidl (on behalf of HERMES Collaboration), 
e-Print Archive: |hep-ex/0405"oT7l talk by N. Makins (on behalf of HERMES Collab- 
oration), Transversity Workshop, Athens, Greece, October 6-7, 2003 

[5] D. Sivers, Phys. Rev. D41 (1990) 83; D43 (1991) 261 

[6] COMPASS Collaboration, P. Pagano, talk delivered at the SPIN2004 Symposium, 
Trieste, Italy, October 10-16, 2004, e-Print Archive: |hep-ex/ 0501035 
COMPASS Collaboration, R. Webb, talk delivered at BARYONS04, Oct 25-29 2004, 



Palaiseau, France, e-Print Archive: hep-ex/0501031 



[7] HERMES Collaboration, M. Diefenthaler, talk delivered at DIS 2005, Madison, Wis- 



consin (USA), April 27 - May 1, e-Print Archive: hep-ex/0507013 

[8] COMPASS Collaboration, V.Yu. Alexakhin et al, Phys. Rev. Lett. 94 (2005) 202002 

[9] J.C. Collins, Phys. Lett. B536 (2002) 43 

[10] M. Anselmino, U. D'Alesio and F. Murgia, Phys. Rev. D67 (2003) 074010 

[11] A.D. Martin, R.G. Roberts, W.J. Stirling and R.S. Thorne, Phys. Lett. B531 (2002) 
216 

[12] S. Kretzer, Phys. Rev. D62 (2000) 054001 

[13] A.V. Efremov, K. Goeke, S. Menzel, A. Metz and P. Schweitzer, Phys. Lett. B612 
(2005) 233 

[14] U. D'Alesio and F. Murgia, Phys. Rev. D70 (2004) 074009 
[15] F. Yuan, Phys. Lett. B575 (2003) 45 

[16] S.J. Brodsky, M. Burkardt and I. Schmidt, Nucl. Phys. B441 (1995) 197 

[17] M. Anselmino, M. Boglione and F. Murgia, Phys. Lett. B362 (1995) 164 

[18] S.J. Brodsky, D.S. Hwang and I. Schmidt, Phys. Lett. B530 (2002) 99 

[19] M. Anselmino, V. Barone, A. Drago and F. Murgia, e-Print Archive: hep-ph/0209073 
A. Drago, Phys. Rev. D71 (2005) 057501 

[20] P.V. Pobilytsa, e-Print Archive: |hep-ph/ 0301236 



11 



[21] A. Bacchetta, A. Schafer and J-J. Yang, Phys. Lett. B578 (2004) 109 

[22] Z. Lu and B.Q. Ma, Nucl. Phys. A741 (2004) 200 

[23] L.P. Gamberg, G.R. Goldstein and K.A. Oganessyan, Phys. Rev. D67 (2003) 071504 

[24] H. Avakian, private communication 

[25] A. Kotzinian, e-Print Archive: hep-ph/0504081 

[26] PAX Collaboration, e-Print Archive: |hep-ex/0505054 

[27] D. de Florian and W. Vogelsang, Phys. Rev. D71 (2005) 114004 

[28] M. Anselmino, M. Boglione, U. D'Alesio E. Leader and F. Murgia, Phys. Rev. D71 
(2005) 014002 

[29] J.C. Collins, Nucl. Phys. B396 (1993) 161 

[30] For a comprehensive recent review paper on transversity, see, e.g. V. Barone, A. 
Drago and P. Ratcliffe, Phys. Rept. 359 (2002) 1 



12 



W 3 
< 



71° 
























7i + (2002-2004) 


HERMES 

not corrected for acce 


: 


ice 


PRELIMINARY 

and smearing 










%■ (2002-2004) 























0.2 0.4 0.6 0.1 0.2 0.3 0.5 1 

z h x B P T (GeV/c) 




Figure 1: HERMES data on ' [7j for scattering off a transversely polarized 

proton target and charged pion production. The curves are the results of our fit. The 
shaded area spans a region corresponding to one-sigma deviation at 90% CL (see text for 
further detail). Predictions for 7r° (upper panel) and kaon (lower panels) asymmetries are 
also shown. 
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Figure 2: COMPASS data 8 on ^ s for scattering off a transversely polarized 

deuteron target and the production of positively (h + ) and negatively (h~) charged hadrons. 
The curves are the results of our fit. The shaded area spans a region corresponding to 
one-sigma deviaton at 90% CL (see text for further detail). 
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Figure 3: The a>dependence of the first moment, according to Eq. (|22jl. of our 
extracted Sivers functions. The solid line is obtained by using the central values of the 
parameters in Table I and the shaded area corresponds to varying the parameters within 
the shaded areas in Figs. 1 and 2. The dot-dashed, dashed and dotted lines show the first 
k± moments of the Sivers functions obtained respectively in Refs. JB], Q3| and |15j . 
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Figure 4: Predictions for A v ^ ~ at COMPASS for scattering off a transversely po- 
larized proton target and the production of positively (h + ) and negatively (h~) charged 
hadrons. The plots in the upper panel have been obtained by performing the integrations 
over the unobserved variables according to the standard COMPASS kinematical cuts PP; 
results with suggested new cuts, Eq. (|27|). are presented in the lower panel. 
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Figure 5: Predictions for * at JLab for the production of 7r + and 7r~ from 

scattering off a transversely polarized proton target. Integrations over the unobserved 
variables have been performed according to the kinematical ranges of Eqs. (|28|) and (|29j) . 
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Figure 6: Predictions for single spin asymmetries in Drell-Yan processes at RHIC, p^ p — > 
£ + £~ X, according to Eq. (j32J) of the text. The lepton pair transverse momentum has been 
integrated in the range < qx < 1 GeV/c; An is plotted as a function of xj? (left) and M 
(right), with integration over the other variable as indicated in the legenda (see text for 
further details). Values of {An} for negative xf are negligible, at RHIC kinematics. 
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Figure 7: Predictions for single spin asymmetries in Drell-Yan processes at GSI, p 1 ' p — > 
£ + £~ X, according to Eq. I|M2J) of the text. The lepton pair transverse momentum has 
been integrated in the range < qx < 1 GeV/c; A jy is plotted as a function of xf (left) 
and M (right), with integration over the other variable as indicated in the legenda (see 
text for further details). The results for the p^ p — > £ + £~ X process are the same. 
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